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The experimental observation of quantum phenomena in mechanical degrees of freedom is dif- 
ficult, as the systems become linear towards low energies and the quantum limit, and thus reside 
in the correspondence limit. Here we investigate how to access quantum phenomena in flexural 
nanomechanical systems which are strongly deflected by a voltage. Near a metastable point, one 
can achieve a significant nonlinearity in the electromechanical potential at the scale of zero point 
energy. The system can then escape from the metastable state via macroscopic quantum tunneling 
(MQT). We consider two model systems suspended atop a voltage gate, namely, a graphene sheet, 
and a carbon nanotube. We find that the experimental demonstration of the phenomenon is cur- 
rently possible but demanding, since the MQT crossover temperatures fall in the millikelvin range. 
A carbon nanotube is suggested as the most promising system. 



I. INTRODUCTION 

The quest towards experimental studies of quantum 
behavior in nanomechanical systems^"- has progressed 
fast in recent years. A breakthrough took place last year, 
when the quantum ground state was demonstrated by the 
Cleland groupP, using a 6 GHz piezoelectric mode in res- 
onance with an electrical quantum system made out of 
Josephson junctions. In general, the difficulty of bring- 
ing mechanical degrees of freedom to the quantum limit is 
due to the challenges posed by several issues. As usual, 
for quantum-mechanical phenomena to become observ- 
able, the thermal energy kbT has to be much lower than 
the characteristic oscillation energy fku. The flexural 
mode frequencies rarely exceed one GHz, and there are 
no such strongly nonlinear mechanical phenomena anal- 
ogous to the Josephson tunneling. Hence, nonlinearity 
which helps to isolate quantum behavior becomes promi- 
nent only by the means of reducing the linear energy, 
with an accompanying reduction in frequency and strin- 
gent requirements for temperature. Also, the mechani- 
cal zero-point vibrations are orders of magnitude smaller 
than the typical length scales encountered in solid-state 
physics. 

Macroscopic quantum tunneling (MQT) refers to 
quantum tunneling in a degree of freedom involving a 
macroscopic number of particles. This is how Joseph- 
son junctions were shown to portray quantum behavior 
more than 20 years ago, namely, by observing the phase 
to escape from the metastable minimum via 
Josephson junctions display strong nonlinearity at the 
zero-point energy scale, and the frequencies typically re- 
side in the range of tens of GHz, and thus the quantum 
limit is encountered at relatively easily attainable tem- 
peratures below one Kelvin. 

A possibility to induce nonlinearity into the mechani- 
cal potential energy was suggested in Refs.lSHIOl In their 
setup, a mechanically induced longitudinal compression 
in a clamped beam would induce a double-well potential, 
and hence a possibility for macroscopic quantum tun- 
neling of displacement of the beam buckling either left 
or right. While measurably high tunneling rates were 



predicted, the compression would have to be adjusted 
with extreme accuracy, and the crossover temperatures 
between thermal activation and MQT would fall in the 
microkelvin regime for beams longer than L ^ 100 nm. 

In the present work, we discuss an alternative way of 
inducing nonlinearity into the etectromechanical poten- 
tial at the scale of zero-point energy. We suggest to 
use a dc voltage to displace a conductive, clamped beam 
or membrane into a deflection close to the critical value 
where it gets pulled in to the electrode. Near this pull-in 
voltage, a metastable minimum appears in the electrome- 
chanical potential. The pull-in event can thus occur via 
a very fundamental physical process, that is, escape by 
quantum tunneling from a metastable state. We see that 
while the control of the deformation requires a high pre- 
cision in the applied gate voltage, the crossover temper- 
ature becomes of the order of mK, and hence is within 
experimental reach. 

In Sec. [n] we first study the behavior of the flexural 
mode eigenfrequency when the beam is influenced by a 
dc electric field from a nearby voltage gate. In Sec. HI 
we calculate in detail the possibility of escaping from the 
metastable minimum via MQT, and consider in particu- 
lar structures based on graphene, or carbon nanotubes. 
Starting from a simplified model, we derive analytical 
expressions for the quantities of interest in the limits of 
very thin or thick beam or sheet. Numerics are used to 
verify the results in a more general setting. 



II. METASTABLE MINIMUM IN THE 
ELECTROMECHANICAL POTENTIAL 

We consider an electromechanical system modeled by a 
doubly clamped beam or sheet, having a length L, width 
W, thickness H, density p and mass m. It is attached 
to the clamps by rigid boundary conditions. An electro- 
static force created by a gate voltage Vg situated at a 
distance d from the beam induces a deformation u{y) to 
it, see Fig. ([I]). Note that we take the direction of the 
thickness H along the direction of deformation. 

Let us write the potential energy of the system. For 
small displacements from the un-stretched position, the 
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to initial tension. Here 
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FIG. 1. (Color online) (a) Schematic picture of the stud- 
ied nanoelectromechanical system. A beam or membrane is 
clamped from opposite ends, and its deformation is controlled 
by an electrostatic force created by a dc gate voltage Vg . (b) A 
metastable minimum in the potential energy near the pull-in 
point is introduced due to the nonlinearities in both electri- 
cal and mechanical energies. The system can escape from the 
minimum via thermal activation (rate Ft), or by macroscopic 
quantum tunneling (rate Fq). 



mechanical energy is quadratic with displacement. Here, 
however, we consider large displacements nearly of the or- 
der d, and take into account the nonlinearity due to the 
elongation of the bearrP^'. This affects the results sub- 
stantially for thin membranes which stretch easily. We 
also include a built-in tension force Tq which exists with- 
out the gate voltage and can be due to fabrication, fn 
addition to the mechanical energy, the total potential en- 
ergy has also an electrical contribution due to the voltage 
bias, given as Vd — —^Cg[u{y)]Vg , where Cg is the ca- 
pacitance between the beam and the gate. The sum of 
all these is 



E. 



'total 




dy 



(1) 

where E is the Young modulus, / the bending moment, 
S the cross section of the beam. We now write the de- 
formation of the beam as u{y) = xox(?/), with xq the 
amplitude and xiu) the mode shape of the deformation. 
Henceforth, we consider the lowest flexural eigenmode, 
whose resonant frequency without the gate voltage is 
u!q — and the mode shape is given by a combi- 

nation of trigonometric and hyberbolic functiond^ (see 
Eq. ([3| for definition of fi) . Scaling the amplitude of the 
deformation with the beam-gate distance, Xq = xg/d, the 
energy becomes 



depend on the shape of the deformation. When using 
the approximatiorP^ for the lowest flexural mode shape 



xiy) 



(1 — cos{2ny)) we get fi = and 



3 



The pull-in phenomenon can be understood as follows. 
The gate capacitance Cg{x), in general, increases towards 
an increasing deflection xq of the membrane. This causes 
an attractive force F = — 9x„Etotai to appear between 
the membrane and the gate. The effective spring con- 
stant 9^Etotai thus contains a positive (three first terms 
in Eq. |2])) and a negative part (last term of Eq. ([2|). 
Therefore it changes sign at a specific voltage Vc depen- 
dent on the geometry, which corresponds to the mem- 
brane getting pulled in to contact with the gate. For 
example, in the simple plate capacitor model, the pull-in 
occurs at a static displacement Xc = \- In an energy 
picture, the minimum of the energy becomes metastable 
when increasing Vg and disappears when Vg = Vc as il- 
lustrated in Fig. 2. 

In order to specify our analysis, we have to choose a 
model for the capacitance. Qualitatively, however, the 
physics of the formation of the metastable minimum re- 
mains model-independent. We consider two cases: (A) 
the parallel plate capacitor, which is that of a membrane 
suspended above a back gate, and (B) a conducting wire 
parallel to a plane. We find carbon-based realizations as 
the most promising, due to their purity, stiffness and low 
mass. The ca ses A and B correspond, respect ively, to 
graphen^i^E^ and a single-wall carbon nanotub d^^ l ^^ l 

First of all, in order to obtain simple analytical esti- 
mates, we assume that the capacitance depends only on 
the amplitude of the deformation xq. This approximation 
is relaxed in the numerical analysis, which takes into ac- 
count the true shape of the deformation (see Appendix) . 
Apart from the thick parallel plate model (section II A 1 1 , 
Xq needs to be computed numerically from the force bal- 
ance equation — VEtotai = 0. 



Etotai =-inujld^xl 



2 j2 

J^d X, 



2 72 2 
-mUJrpd Xn 



1/2 
•^5 



Cg[u{y)]dy, 



(2) 



A. Parallel plate model 



Using the parallel plate capacitance model Cg 
d(i^o) ' energy becomes 
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with an increasing voltage, whereas the fourth term, elec- 
trical nonlinearity, has the o pposi te effect. In exp eri- 
ments with carbon nanotube d^^ * ^^ ! or graphene^^ES, the 
frequency has typically been observed to go up, but for 
microfabricated metallic resonators, the frequency has 
decreasecP^^. We now study the behavior of the fre- 
quency when we get close to the pull-in for two different 
regimes. In the first regime, we consider the case where 
the bending energy is large compared to the stress en- 
ergy, i.e^>l 4^ jj = f3 <^ 1. This is the case when 
the variation of the frequency coming from the induced 
deformation is small, e.g., when the beam is close to the 
gate or when the beam is thick. We then consider the 
opposite regime where it is the first term of Eq. Q which 
can be neglected, e.g., for a thin beam or a beam far from 
the gate {j^ > 1). 




1. Thick membrane, jS — ^ 1 

For a beam close to the gate, the electrical potential is 
dominant, and the third term on the rhs of Eq. Q can 
be ignored. This also holds if the built-in tension is large, 
since a larger voltage is then needed in order to obtain a 
certain xq. Scaling the energy with the bending energy 
we define Ethick = Etotai/(^TOWod^): 



FIG. 2. (Color online) Schematic picture of the pull-in phe- 
nomenon in a suspended membrane or beam: when increas- 
ing the dc voltage, the amplitude of the deformation increases 
until the beam pulls into contact with the gate at a certain 
critical voltage Vc. 
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Expanding Eq. ([sj up to the third order for small varia- 
tions X around the equilibrium amplitude xoiVg) we ob- 
tain 
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with 
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An inspection of Eq. ([t]) reveals how the frequency de- 
pends on the gate voltage. The third term on the rhs, the 
mechanical nonlinearity, tends to increase the frequency 



Ethick 



V (1- 2^0)4 ■ 
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with y2 ^ 



eoWL 



The critical displacement at which 



pull-in point occurs is given by Xc = l/S,!^ and the 



critical voltage is Vc/V 
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Using the 



plate-capacitor model, the position-dependent frequency 
and the third-order coefficient near the pull-in point are 
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For example, a graphite resonator with p = 2 g/cm'^, 
i = ly = 0.5 Aim, = 100 nm, d = 10 nm and = 30 
GPa, would have K = 25.9 V and a;o/27r = 1.6 GHz. 
A numerical calculation of the behavior of the frequency 
with respect to the gate voltage is shown in Fig. 3. In 
this figure and Figs. 4-7, d has been chosen so that the 
accumulated strain in the pull-in position is at most half 
of the tensile strength.l^ 



2. Thin membrane, /3 = j| S> 1 

The electrical potential term in Eq. ([ij is insignificant 
at small deflections, if the beam is thin or far from the 



4 



(a) 



(a) 



8 O 

3l 3 



0.8 



0.6 



0.4 



0.2 



1 r 




-Numerical calculation 
-Analytic calculation 



_l I I L_ 



0.1 0.2 0.3 0.4 0.5 0.6 

V /V 



(b) 

0.04 

(M h|(mo0.03 
3l 3 

0.02 
0.01 




"Numerical result 



1.5 1 0.5 



FIG. 3. (Color online) Behavior of the first flexural mode fre- 
quency with respect to the applied gate voltage scaled with 

V — \J '^^'^'^ for a thick membrane having To = Q, L = W 
and /3 = 0.85. (a) Black solid curve: analytical calculation 
supposing a constant mode shape, red dashed curve: full nu- 
merical calculation with the exact mode shape, (b) Behavior 
of the frequency close to the pull-in point. 
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FIG. 4. (Color online) As Fig. [s] but for a thin membrane 

with {V'f = ^^1^, L = W and /3 = 283. (a) Full numerical 
calculation using the exact mode shape, with two different 
built-in tensions. The blue dashed curve: the initial tension 
is To = while for the red solid curve it is To = 2 ■ lO'^ f| . (b) 
Fit of the numerical result close to the pull-in point with the 
equation a/3^V3 [l - {Vg/Vcf]"". The fit gives a = 0.7 and 
m = 0.5 which confirm the power law in Eq. ( 13 1. 



gate. The beam then becomes stiffer and the frequency 
goes up with the gate voltage as shown in Fig 4. 

At larger deflections, however, the electrical term be- 
comes dominant and the metastable potential minimum 
becomes shallow and eventually disappears. In this 
case we may neglect the small constants uiq and in 
Eq. The pull-in point is obtained when the elec- 

trical and mechanical terms in Eq. ([t]) cancel one an- 
other. Scaling the energy with the stress energy we define 
Ethin = Etotai/^mw^d^): 



^thii 



1 



12x^ - 2 



V,V 1 



4xr ' 



V*J (l-.To)3 
2 



1 



V*J {l-xoY 



(12) 



with (V* 



enWL 



We obtain the critical displace- 



ment Xc — 3/5, and the critical voltage^^ Vc/V* = 



II w |. Near the pull-in point, we obtain from Eqs. (7 



8) 



a 



54 
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2 fJL 
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Contrary to a thick membrane, Eq. (Ill, the nonlinearity 
a depends on jP . Therefore the nonlinearity is linked to 
the elongation of the beam. For example, taking a single- 
layer graphene sheet with p — 2 g/cm^, L = = 0.5 
/im, H = 0.3 nm, d = 85 nm and E = 1 TPa, we get 
Vc = 17.8 V and wo/27r = 27.6 MHz. 
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B. Wire parallel to a plate 

The parallel plate capacitor model is a good approxi- 
mation for membrane resonators, but another interesting 
case to consider is a suspended beam, such as a thin car- 
bon nanotube. With the capacitance Cg — inprf/H(\^^o)] ' 
H being in this case the radius of the beam, the energy 
of the system is 



As in the previous section we make an expansion of 
Eq. ( 14 ) up to the third order for small variation x around 



the equilibrium amplitude xo(Vg) and obtain 
1 



E. 



with 
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a = 2muj'^d^xo 

Vfvreo (ln2(2(l-a:o)/3) 



l)2 1n3[2(l-a;o)/3] 

31n(2(l-a:o)/3) 



(16) 
(17) 



3(a;o - l)Mn*(2(l - .To)/3) 

The behavior of the frequency is similar to the one 
found in the previous section. Below we only consider 
the case of a thin tube ji ^ ^ which best corresponds to 
carbon nanotube^^H 



1. Thin tube, /3 = |. > 1 

As in the previous section, we may neglect the small 
constant Wq term in Eq. (16). The pull-in point is 



obtained when the electrical and mechanical terms in 



Eq. ( 16 1 cancel one another. Neglecting in the derivatives 



the xq dependence in the logarithms, we obtain the criti- 
cal displacement Xc = 2(^4-2\n2B\ ^'^'^ critical voltage 



2(1+2 In 2^) 
^2 _ 27/3"ln2/3°(2+ln2/?) ^j^j^ y2 
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pull-in point, we obtain from Eqs. ( 16|17 1: 
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Near the 



(18) 



-C, 



where 



C 



In(^) 



321n^(/3) 



Similarly to the thin membrane, the nonlincarity is 
due to the elongation of the beam, but there is a differ- 
ence in the behavior of the frequency close to the pull-in. 



Here, decreases linearly with the voltage, instead of 



being proportional to a square root of Vg as in Eq. ( 10 1 



Therefore, the pull-in for a beam is expected to occur in 
a "smoother" way than for a membrane. The results for 
the frequency are summarized in Fig 5. Although the 
qualitative behavior of the frequency is mostly captured 
by the analytical calculation, very close to the pull-in we 
again capture the square root behavior (Fig. 5b). For 
example, with a single-wall carbon nanotube resonator 
with p = 2 g/cm^, L = 0.5 /^m, H = 1 nm, d = 85 nm 
and S = 1 TPa, we have = 6.2 V and Wq = 23 MHz. 



(a) 




FIG. 5. (Color online) The fundamental flexural mode fre- 
quency with respect to the applied gate voltage with V'^ = 



for a single-wall carbon nanotube with /? — 85. (a) 



2 ,3 
maj^ a 

The blue dashed curve is without an initial tension (To = 0), 
while the red solid curve is with To = 3.1 • 10^f|. (b) Fit 
of the numerical result with the equation a/3^(l — (Vg/Vc)'^)"^ 
indicates a power law with a = 0.32 and m = 0.5. 



6 



III. PROSPECTS OF OBSERVING MQT 

At a temperature T, the rate of thermally activated es- 
cape from the metastable minimum in Fig. 3 is given by 
the usual Arrhenius law Ft = ciJa;/(27r) exp (— AT^/fc^T) 
where w^, is the frequency of the resonator. The tun- 
neling rate from the ground state is given bj^^H Fq = 
loq exp (— 1^ ^) , where N = is the number of states 

in the potential with height AV, and loq = 6uJx\/GN/tt. 
The crossover temperature Tq below which quantum tun- 
neling from the ground state dominates over thermal es- 
cape is given by equating the Arrhenius law with the 
quantum tunneling rate, 



10' 



To 



5 hujx 



1 



36 kh 



1 



i^in 12 



(19) 



In addition to giving the crossover between the ground- 
state tunneling and thermal activation, this is also a good 
criterion for the system to be in the quantum limit in 
the sense that thermal population becomes negligible. 
Another important quantity is the number N of bound 
states in the metastable minimum. From Eq. ^ we ob- 
tain 



N 



AV 

hlOrr 



1 d^m^uj^ 
54 ha'^ 



(20) 



In order to maximize the quantum tunneling rate we need 
to have ^ 1. For a membrane in the case where /3 <C 1 
the frequency at which the number of states is close to 1 
satisfies 



1^ 

UJQ 



3^ 1 



while when /3 ^ 1 it is 



UJQ 



1100 



(21) 



(22) 



with A^o = mLOocP/fi- From these equations we can see 
that the mass of the resonator has to be the lowest pos- 
sible in order to maximize the frequency, which would 
then lead to the highest crossover temperature, and hence 
the best experimental prospect of observing MQT. Also 
Eq. (22) shows that a system having a large /3 is favor- 



able since it yields a higher frequency when the num- 
ber of states is close to 1. The most suitable materials 
from both of these points of view, having a low mass and 
possibility for a high ratio, are graphene and carbon 
nanotubes. In Figs.|6]and [7]we plot a full numerical cal- 
culation of the crossover temperature and the MQT rate 
for graphene and CNT resonators, which show the con- 
flicting requirements of attaining a high frequency and 
simultaneously having a measurable tunneling rate. 

MQT could be verified by measuring the reduction 
of the escape rate as a function of temperature, and 
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FIG. 6. (Color online) Macroscopic quantum tunneling 
(MQT) in a thin membrane, /3 2> 1. Right hand axis in 
green: the quantum-classical crossover temperature, plotted 
as a function of departure of the dc gate voltage from the crit- 

ical voltage with V — , 4 . Left hand axis in blue: the 
rate of MQT of the displacement out of the metastable min- 
imum. The plots are for a representative sample consisting 
of a single-layer graphene sheet with H = 0.3 nm, L = 5 fim, 
W = 2//m, E = 1 TPa, p = 2 g/cm^, and d = 85 nm. The 
lowest number of states is A ~ 2. 




(Vc - Vg)/V 



FIG. 7. (Color onHne) As Fig. [g] but for a carbon nanotube, 
with H ^ 1 nm, L = 0.5 //m, E = 10^ GPa, p = 2 g/cm^, and 

d = 85 nm. Here V^ = T^. 



observing its saturation at the crossover temperatures 
predicted above. An escape event would easily be de- 
tected as a large change of capacitance once the mem- 
brane g ets pu lled in, for instance, by using the dispersive 
methodJi^ElHlll which do not otherwise excessively dis- 
turb the system. In order to repeat the experiment, the 
gate voltage would be reset to zero, and another gate on 
the opposite side could be used to pull the membrane 
from the Van der Waals attraction. 

A possible experimental verification is challenged by 
the fact that reaching a small number of states requires 
a high precision on the applied gate voltage, and a small 
change in the applied gate voltage leads to a large varia- 
tion in the frequency in this regime. The wire parallel to 
a plate model shows that the variation close to the pull- 
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in point is smoother than for the parallel plate model, 
and in this respect, carbon nanotubes seem to be good 
candidates for observing MQT. 



IV. CONCLUSIONS 

Macroscopic quantum tunneling (MQT) is a fascinat- 
ing topic which allows one to investigate the quantum- 
classical transition regime, where frequencies of collective 
degrees of freedom become comparable to temperature. 
We have introduced a model system of a mechanical de- 
gree of freedom trapped into a metastable state formed 
by a conductive beam or membrane suspended on a volt- 
age biased back gate. Via MQT, the fictitious particle 
can escape from the metastable state and thereby be 
pulled into contact with the gate. The setup might serve 
as a means of observing mechanical MQT, in a fashion 
analogous to how quaiitum behavior in Josephson sys- 
tems was first observed^-"^. 

Here we discuss the possibility of observing MQT in 
nanomechanical resonators for two capacitance models, 
those of a parallel plate capacitor, and the wire paral- 
lel to a plate, corresponding, for example, to suspended 
graphene or carbon nanotube. One should have the num- 
ber of states in the metastable minimum close to one in 
order to obtain a sufficient quantum tunneling rate, and 
at the same time, maintain a high frequency to maxi- 
mize the crossover temperature. The highest crossover 
temperatures are obtained in a system with low mass, 
low density, high Young's modulus, and a high ratio 
d/H which allows to increase the frequency by increased 
strain. These criteria point towards carbon-based real- 
izations. 

We conclude that while the predicted crossover tem- 
peratures in the mK-range are several orders of mag- 
nitude higher than for the buckled beam studied pre- 
viously, they are still experimentally demanding, and 
barely within reach of standard dilution refrigerator tech- 
niques. However, one may use electrical cooling tech- 
niques where the nanoresonator is coupled to higher- 
frequency electrical resonator in order to cool the lowest 
mode in questio n dow n to temperatures much lower than 
the environmentP^'29l_ ^j^^g foresee the experimental 
verification challenging, but possible in the future. 
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APPENDIX 

For the numerical calculations, we assume that there 
is a small deviation Su{y,t) = Mi(j/)e"^* from the static 
deformation uo{y) such that the total deformation can 
be written in the form 



u{y,t) = uo{y) + 6u{y,t). 



(23) 



Introducing these expressions in the Euler-BernouUi 
equation, neglecting the terms which are O(w^) and tak- 
ing the parallel plate capacitance model leads to the 
equation for DC defiection uq 



dy^ 



dy'^ 
1 



(l-"o)^ 



0. (24) 



The AC part ui{y,t) satisfies an eigenvalue equation for 

Ml 
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Here, 



= —Tn 



/3 = 



EI 
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H 
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EI 



(25) 

(26) 
(27) 

(28) 

(29) 



Writing uq = YJ- a^Xiiu) with Xiiu) the i*-^ flexural 
eigenmode,'^^ we solve ( |24| ) using the Galerkin method, 
rewriting ( 25 ) in the eigenmode space (xi-space) and use 
the solution found for uq to compute uj and ui. This 
first step allows us to find the behavior of the frequency 
with respect to the applied gate voltage. This is shown 
in Figs. [3] and [4] in the thick and thin membrane regime 
limits, for /3 < 1 and /3 > 1, respectively. We then use 
Mo(j/) and ui{y,t) and use Eq. ([ij to plot e{uQ + xqUi) 
for a particular gate voltage, xq being the amplitude of 
the AC deformation. This step allows us to compute 
the height of the metastable potential AV. From these 
results it is possible to compute the number of states and 
thus the quantum tunneling rate as discussed in Sec. HI 
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